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PREFACE 


This  Memorandum  presents  an  algorithm  for  maximizing  the  flow  from 
one  node  to  another  through  a  network  whose  arcs  and  nodes  have  limited 
capacities.  It  Is  a  product  of  RAND  and  Air  Force  Interest  In  network 
flows  and  Is  currently  being  used  In  an  Interdiction  study  which  Is  the 
subject  of  a  forthcoming  RAND  Memorandum.  If  the  reader  is  Interested 
in  an  extensive  description  of  the  foundations  of  network  theory,  he 
should  read  L.  R.  Ford  and  D.  R.  Fulkerson,  Flows  In  Networks.  The  RAND 
Corporation,  R-375-PR,  December  1960. 

This  study  should  be  useful  to  persons  in  the  Operations  Analysis 
Office,  Hq.  USAF  (AFGOA) ,  the  Intratheater  Transportation  Study  Group, 
and  others  who  deal  with  Interdiction  anl  transportation  problems. 
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SUMMARY 

In  many  actual  network  flow  situations,  nodes  as  well  as  arcs  have 
limited  capacities.  This  Memorandum  presents  for  such  a  network,  an 
algorithm  for  maximizing  flow  from  a  source  node  to  a  sink  node.  The 
algorithm  allows  us  to  treat  these  situations  without  introducing 
artificial  arcs  and  nodes,  as  has  been  done  in  the  past.  Eliminating 
the  artificial  arcs  and  nodes  simplifies  network  analysis  since  it 
always  results  in  half  as  many  nodes,  as  well  as  less  than  half  as 
many  arcs  if  the  original  arcs  are  undirected. 

In  addition,  the  following  generalization  of  Ford  and  Fulkerson's 
max-flow,  min-cut  theorem  is  presented  and  proven.  (See  L.  R.  Ford 
and  D.  R.  Fulkerson,  Flows  in  Networks.  The  RAND  Corporation,  R-375- 
PR,  December  1960.)  Consider  two  subsets  of  the  nodes,  X  and  Y,  whose 
union  is  the  set  of  all  network  nodes  and  such  that  the  source  node  is 
a  member  of  X,  and  the  sink  node  is  a  member  of  Y.  Then,  forming  a 
cut  set  separating  the  source  and  sink  are  the  nodes  in  the  intersection 
of  X  and  Y,  and  the  set  of  all  arcs  (i,  j) ,  such  that  t  is  a  member 
of  X  -  Y,  and  j  is  a  member  of  Y  -  X.  Letting  a  cut  set's  value  be 
the  sum  of  the  capacities  of  all  its  arcs  and  nodes,  it  follows  that 
the  maximum  flow  is  equal  to  the  minimum  value  of  all  cut  sets  separating 


the  source  and  sink. 
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I.  INTRODUCTION 

The  problem  dealt  with  in  this  Memorandum  is  one  of  maximizing 
the  flow  through  a  network  in  which  the  nodes  and  arcs  have  limited 
capacity.  Specifically,  a  network  consisting  of  arcs  and  nodes  is 
given.  The  nodes  are  points  where  two  or  more  arcs  meet,  and  the 
arcs  are  line  segments  connecting  two  nodes.  Each  arc  has  a  positive 
capacity  representing  the  maximum  amount  of  flow  that  may  pass  across 
it,  and  each  node  has  a  positive  capacity  representing  the  maximum 
amount  of  flow  that  may  enter  or  leave  it.  If  the  network  is  directed, 
each  arc  also  has  a  direction  indicating  the  direction  along  which  flow 
is  allowed  to  pass.  For  undirected  networks,  flow  may  pass  in  either 
direction  on  an  arc.  Figure  1,  page  b ,  shows  an  example  of  a  directed 
network  where  it  is  desired  to  find  the  maximum  amount  that  may  past 
from  a  node  known  as  the  source  to  one  called  the  sink,  subject  to  the 
arc  and  node  capacity  constraints. 

A  previous  study  of  this  problem  split  each  node  into  two  nodes 
* 

and  an  arc.  In  the  case  of  undirected  networks,  this  method  also 
requires  replacing  each  arc  by  two  directed  arcs.  The  algorithm  pre¬ 
sented  here  should  allow  substantial  savings  in  computer  time  and  core 
storage  for  undirected  networks  and  may  also  allow  good  savings  for 
directed  networks. 

Section  II  describes  the  linear  program  that  maximizes  flow  through 
the  network.  Section  III  illustrates  how  cut  sets  are  used  to  solve  the 
problem,  and  Sec.  JV  gives  the  algorithm  for  finding  a  maximum  flow 

L.  R.  Ford  and  D.  R.  Fulkerson,  Flows  in  Networks,  The  RAND 
Corporation,  R-375-PR,  December  1960. 
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pattern.  Section  V  contains  proofs  that  the  linear  program's  constraint 
(conservation  at  the  nodes,  arc  capacities,  and  node  capacities)  hold. 
Section  VI  shows  how  the  maximum  flow  algorithm  terminates  and  that  the 
maximum  flow  equals  the  minimum  cut,  while  Sec.  VII  demonstrates  the 


algorithm  with  an  example. 
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II.  LINEAR  PROGRAM  FORUMULATION 

The  problem  of  maximizing  flow,  v,  through  a  network  with  node 
and  arc  capacities  can  be  formulated  as  a  linear  program,  tut  c(i) 
be  the  capacity  of  node  l;  c(i,  j)  the  capacity  of  arc  (i,  J)  ;  x(i,  j) 
the  flow  on  arc  (i,  j) ;  A(i)  the  set  of  nodes  to  which  arcs  are  directed 
from  node  i;  B(i)  the  set  of  nodes  from  which  arcs  are  directed  to  i;  S 
the  source  node;  and  S  the  sink.  The  linear  programming  formulation  of 
the  network  problem  is  as  follows. 

Maximize  v  subject  to: 

v,  x(i,j)  S  0  all  i  , j 

I  x(i,j)  *  c(i,j)  all  <i,J) 

II  a)  £  x(i,j)  -s  c(i)  i^S 
jeA(i) 

b)  £  x(j,S)  <  c(S) 
jcB(S) 

Ilia)  £  x(s»j)  " 

jcA(S)  jcB(S) 

b)  £  x(i  ,  J)  -  £  x(j,i)  =  0  iilS.S 

jcA(i)  jeB(i) 

c)  £  x<S,j)  -  £  x(j,5)  -  -v  . 

j«A(S)  j  eB(S) 

Relationship  I  represents  the  arc  capacities  and  II  the  node 
capacities.  Note  that  for  all  nodes  other  than  the  source  and  sink. 

For  undirected  networks,  c(i,  j)  ■  c(J,  i)  for  all  (i,  J) . 


£  x(j  ,S)  =  v 
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the  nods  capacity  may  be  expressed  either  in  terms  of  the  flow  into 
or  the  flow  out  of  the  node  since  these  two  quantities  are  equal*  For 
the  source  and  sink,  these  two  quantities  are  unequal  and  the  node  ca¬ 
pacity  must  be  expressed  in  terms  of  the  greatest  of  these  two  quanti¬ 
ties*  Thus  the  source  node  capacity  limits  flow  out  of  the  source  while 
the  sink  capacity  limits  flow  into  the  sink*  Relationships  Ill  are  bal¬ 
ance  equations  showing  that  the  net  flow  out  of  the  source  and  the  net 
flow  into  the  sink  both  equal  the  value  of  the  flow  and  that  the  net 
flow  into  or  out  of  any  other  node  is  zero*  Note  that  any  one  of  the 
balance  equations  may  be  dropped  since  they  sum  up  identically  to  zero* 
Physically,  this  dependency  expresses  the  fact  that  every  unit  of  arc 
flow  enters  a  node  and  leaves  a  node.  For  the  network  of  Fig.  1,  the 
linear  program  would  be,  maximize  v  subject  to: 

v,  x(i, j)  s  0  all  i, j 
x(S,a)  £  5 
x(S,b)  s  1 
x(S,c)  *  6 
x(a,b)  £  3 
x(a,S)  £  5 
x(b,c)  £  4 
x(b,S)  £  7 
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x(c ,  S)  ^  2 

x(S  ,  a)  +  x(S,  b)  +  x(S ,  c)  £  12 
x(a ,  b)  +  x(a,  S)  s  10 
x(b ,  c)  +  x(b  ,  S)  ^  8 
x(c ,  S)  <■  7 

x(a,  S)  +  x(b ,  S)  +  x(c  ,  S)  <  11 

x(S ,  a)  +  x(S ,  b)  +  x(S  ,  c)  -  v  =  0 

x(a,  b)  +x(a,  S)  -  x(S,  a)  =  0 

x(b,  c)  +  x(b,  S)  -  x(S,  b)  -  x(a,  b)  =  0 

x(c,  S)  -  x(S,  c)  -  x(b,  c)  =  0 

-x(a,  S)  -  x(b  ,  S)  -  x(c  ,  S)  +  v  =  0  # 

A  set  of  nodes  a. . a  is  called  a  chain  from  a  to  a 

in  in 

if  either  (at ,  ai+1)  or  (ai+1>  at)  is  an  arc  for  1  a  1 . .  -  1. 

If,  in  this  chain,  (a^,  a^)  is  an  arc  for  eacl'  ^  the  chain  is  a 
directed  one.  A  chain  in  which  all  the  a^  are  distinct  is  a  path. 
A  chain  in  which  a^  ■  aQ  is  called  a  cycle.  Also  paths  and  cycles 
that  come  from  directed  chains  are  directed  paths  and  cycles. 
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Ford  and  Fulkerson  show  that  any  set  of  x(i,  j)  a  0  which  satisfies 
the  balance  equations  may  be  expressed  as  the  sum  of  flows  along  directed 
paths  from  S  to  S  and  directed  cycles  where  the  sum  of  flows  along  the 
paths  is  v.  Furthermore,  if  both  arcs  (i,  J)  and  (J,  i)  exist,  any  fea¬ 
sible  flow  value  v  may  be  obtained  with  non-zero  flow  on,  at  most,  one 
of  these  arcs,  by  replacing  x(i,  j)  with  max  [x(i,  j)  -  x(J ,  i) ,  0], 

* 

Ibid. 
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m.  CUT  SETS 

The  solution  method  presented  in  this  Memorandum  will  take  advan¬ 
tage  of  the  properties  of  a  cut  set,  defined  below.  This  definition  is 
a  generalization  of  the  cut  set  definition  given  in  Ref.  1  for  solving 
network  problems  with  arc  capacities. 

Definition:  Given  a  set  of  nodes  N  and  arcs  A,  consider  two  sub¬ 
sets,  X  and  Y,  of  the  nodes  such  that  X  U  1  ■  N,  If  S  e  X  and  S  e  Y, 
then  the  set  of  nodes  in  X  0  Y  together  with  the  set  of  arcs  (i,  J)  , 
such  that  1  c  X  -  Y  and  j  e  Y  -  X,  form  a  cut  set  separating  S  and  S. 

The  value  V(Xt  Y)  of  this  cut  set  is  equal  to  the  total  capacity  of  its 
arcs  and  nodes. 

If  S  and  S  are  the  source  and  sink,  respectively,  it  follows  that 
since  the  source  is  in  X  and  the  sink  is  in  Y,  any  path  from  source  to 
sink  must  use  at  least  one  arc  or  one  node  in  the  cut  set.  Thus  the 
maximum  flow  cannot  exceed  the  value  of  any  cut  set.  In  particular, 
it  cannot  exceed  the  minimum  value  of  all  cut  sets.  It  will  be  shown 
later  on  that,  as  in  networks  with  arc  capacities  only,  the  maximum  flow 
equals  the  minimum  cut.  Some  pertinent  relationships  between  flows  and 
cuts  follow. 

Lemma  1.  Suppose  a  flow  pattern  con.*  isting  of  one  unit  of  flow 
along  a  directed  chain  from  the  source  node  S  to  node  a  is  given.  Sup¬ 
pose  also  that  (X,  Y)  is  a  cut  set  and  S  e  X  -  Y.  Then  for  i  e  X  -  Y, 

J  e  Y  (note  that  X  -  Y  and  Y  are  complements)  : 

5>U>  j)  -  Ex(i.  O-l.  a  «  Y 

0,  a  e  X  -  Y  • 

Proof.  A  quick  check  shows  that  the  theorem  holds  if  the  chain 
consists  of  one  arc.  Suppose  it  holds  for  k  arc  chains.  Consider  an 
arbitrary  k  +  1  arc  chain,  S,  a^  ...  ,  a^+1*  There  are  four  possible 
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cases  Co  consider. 

Case  l!  «k,  ak+l  a  Y. 

Since  by  hypothesis  the  theorem  holds  for  k  arc  paths, 
the  above  quantity  must  be  one  if  one  neglects  the  unit 
of  flow  on  arc  (a^,  a^+j)  .  Furthermore,  the  unit  of  flow 
on  this  arc  contributes  nothing  to  the  quantity  in  ques¬ 
tion,  so  it  must  still  be  one. 

Case  2:  a^  e  Y,  a^+^  c  X  -  Y. 

Neglecting  the  flow  on  the  last  arc  yields  a  value  of  one 
for  the  quantity,  and  the  unit  of  flow  on  arc  (a^,  a^j) 
reduces  it  to  zero. 

Case  3:  afc  e  X  -  Y,  ak+1  e  Y. 

Neglecting  the  flow  on  the  last  arc  yields  zero.  Adding 
the  flow  on  the  last  arc  adds  one  to  the  quantity  bringing 
it  to  one. 

Case  4;  a^,  afc+1  e  X  -  Y. 

Neglecting  the  last  arc  yields  zero  for  the  quantity. 
Adding  the  unit  of  flow  on  the  last  arc  does  not  change 
this. 

Since  the  theorem  is  preserved  in  all  four  cases,  the  theorem  itself 
follows  from  induction. 

Corollary  1.  Suppose  a  flow  pattern  is  given  consisting  of  one 
unit  of  flow  along  a  directed  cycle.  Then  the  quantity  referred  to  in 
Lemma  1  is  zero. 

Proof.  Suppose  the  cycle  contains  a  node  in  X  -  Y.  Designating 
one  such  node  as  S,  one  has  a  directed  chain  from  S  to  S  for  which  the 
quantity  in  question  must  be  zero  by  Lemma  l.  On  the  other  hand,  if  all 
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nodes  of  the  cylce  are  in  Y,  then  all  x(i(  j)  in  the  summations  refer¬ 
red  to  are  zero  and  the  theorem  still  holds. 

Theorem  1:  Suppose  a  flow  pattern  is  given  satisfying  the  balance 
equations  (III)  of  Sec.  II.  Suppose  also  that  S  e  X  -  Y  and  S  e  Y. 

Then 

£*(i,  J)  -  2>(J,  t)  -  v,  i  e  X  -  Y,  j  e  Y  . 

Proof:  Multiplication  by  v  and  Lemma  1  show  that  the  theorem 
holds  if  the  flow  pattern  consists  of  v  units  of  flow  along  a  directed 
chain  from  S  to  S.  Similarly,  Corollary  1  shows  it  holds  if  the  pattern 
consists  of  units  of  flow  only  along  one  directed  cycle.  The  theorem 
Itself  follows  since  the  flow  pattern  can  be  decomposed  into  units  of 
flow  along  directed  chains  from  S  to  S  and  directed  cycles  such  that  the 
sum  of  these  units  on  the  chains  is  v. 

Theo rem  2 ;  The  maximum  flow  cannot  exceed  the  value  of  any  cut 
set  separating  S  and  S. 

Proof :  Suppose  S  e  X  -  Y.  By  Theorem  l,  v  ^  x(i,  J),  1  e  X  •  Y, 
j  e  Y.  However,  the  right  hand  side  of  the  above  inequality  is  equal  to 
£x(i,  j)  +  £  x(i,  k)  ,  i  e  X  -  Y,  j  e  Y  -  X,  k  e  X  n  Y. 

The  first  summation  term  in  this  expression  cannot  exceed  the  total  arc 

capacities  in  (X,  Y)  ,  and  the  second  cannot  exceed  the  total  node  capac¬ 
ities  in  (X,  Y)  .  Thus  the  theorem  holds  for  S  e  X  -  Y.  If  S  e  X  fi  Y, 

(X,  Y)  contait.s  node  S  and  therefore  it  suffices  to  show  that  v  £  c(S) . 

To  see  this,  attach  an  artificial  node  S'  and  an  artificial  arc  (S',  S)  , 
each  with  infinite  capacity.  Let  S'  be  the  source  and  S  the  sink  of  the 
new  network.  Note  that  any  feasible  flow  (one  that  satisfies  the  con¬ 
straints  of  the  linear  program  in  Sec.  II)  can  be  converted  into  a 
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feasib  w  for  the  new  network  by  setting  x(S,  S')  .  v.  Thus  the 

maximum  now  of  the  orlgi  il  network  cannot  exceed  the  maximum  flow 
of  the  new  one.  Furthermore,  letting  X  consist  of  the  nodes  S'  and  S, 
and  Y  consist  of  all  nodes  except  S',  gives  a  cut  set  whose  value  is  at 
least  as  great  as  the  maximum  flow  of  the  network  (since  S'  e  X  -  Y)  and 
is  equal  to  the  capacity  of  node  S. 

These  relationships  will  be  used  to  show  that  the  algorithm  pre¬ 
sented  in  Sec.  IV  yields  a  maximum  flow. 
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IV.  A  MAXIMUM  FLOW  ALGORITHM 


The  algorithm  for  finding  a  maximum  flow  pattern  starts  out  with 
a  flow  of  zero  on  all  arcs.  Then  the  algorithm  attempts  to  find  a 
chain  from  source  to  sink  wlr.h  unused  capacity  on  ail  of  its  arcs  and 
nodes.  If  one  is  found,  as  much  flow  as  possible  is  sent  across  it. 

Then  an  attempt  is  made  to  find  a  new  chain,  until  finally,  none  can 
be  found  and  the  algorithm  terminates.  The  chains  found  will  not  neces¬ 
sarily  be  directed.  This  means  that,  in  general,  sending  flow  along 
the  chain  increases  the  flow  on  some  arcs  and  decreases  it  on  others. 

The  chains  found  may  contain  cycles.  It  will  be  shown,  however,  that 
while  nodes  may  appear  more  than  once  in  a  chain,  they  may  not  appear 
more  than  twice,  and  that  arcs  may  appear  only  once.  Furthermore,  only 
a  saturated  node  may  appear  twice  with  the  included  cycle  reducing  the 
node  flow  to  balance  out  in  Increase  by  the  rest  of  the  chain.  Flow 
can  be  pushed  across  a  chain  until  either  the  flow  on  one  of  its  arcs 
or  nodes  reaches  capacity  or  until  the  flow  on  one  of  its  arcs  is  reduced 
to  zero. 

The  attempt  to  find  a  chain  from  source  to  sink  is  a  labeling  tech¬ 
nique  in  which  the  presence  of  a  label  at  a  i  ode  indicates  that  a  chain 
with  unused  capacity  has  been  found  from  the  source  to  it.  Initially 
the  source  is  labeled  and  all  other  nodes  arc  unlabeled.  Attempts  are 
made  to  label  the  unlabeled  nodes  by  examining  the  labeled  ones.  If 
the  sink  becomes  labeled,  units  of  flow  are  sent  along  the  source-sink 
chain  found  and  the  labeling  process  is  repeated  in  an  attempt  to  find 
a  new  source-sink  chain.  The  algorithm  terminates  when  no  such  chain 


can  be  found. 
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The  label  at  an  arbitrary  labeled  node,  a,  wLU  be  designated  by 
the  three-component  vector  L(a)  *  (a^,  a2*  a3^  ’  t^ie  Individual 

components  being  referred  to  as  L^(a)  ,  L^(a)  ,  and  L^(a).  In  addition, 
each  node,  a,  labeled  or  unlabeled  will  have  a  quantity  c(a)  which 
represents  Its  unused  capacity.  The  meaning  of  the  label  components 
are  as  follows:  If  L^(a)  =  b+,  then  a  chain  to  node  a  has  been  found 
along  which  there  are  at  least  L^(a)  units  of  unused  capacity  except 
possibly  at  node  a,  and  the  last  arc  along  this  chain  is  (b,  a).  If 
L^(a)  =  b  then,  as  before,  a  chain  has  been  found  from  the  source  to 
node  a  along  which  there  are  at  least  L^(a)  units  of  unused  capacity; 
however,  in  this  case,  the  last  arc  is  (a,  b)  .  Note  that  if  units  of 
flow  were  sent  along  the  chain,  an  increase  in  the  quantity  x(b,  a) 
would  result  in  the  former  case,  while  a  decrease  in  the  quantity  x(a,  b) 

would  result  in  the  latter  case.  The  quantity  L^a)  ,  if  it  exists,  can 

only  take  the  form  b  and  has  the  same  meaning  as  L^(a).  It  will  be 
assigned  only  if  no  units  may  be  sent  to  a  along  the  chain  designated 
by  L^(a)  without  violating  the  node  capacity  at  node  a.  This,  of  course, 
can  only  happen  if  L^(a)  is  of  the  form  b+ ,  and  c(a)  =  0.  If  L^(a)  does 
not  exist,  it  will  be  designated  by  the  symbol  *.  When  L^Ca)  4  *,  in 
which  case  two  such  chains  to  a  have  been  found,  the  quantity  L^(a) 
will  be  such  that  L^(a)  units  may  be  sent  along  either  chain. 

If  the  above  conditions  hold,  note  that  if  i  is  labeled  and  j  is 
unlabeled  with  x(j  ,  i)  >  0,  then  j  may  be  labeled,  since  any  chain  to 

node  i  with  unused  capacity  plus  a  flow  reduction  on  arc  (j ,  i)  is  a 

chain  to  j  with  unused  capacity  (Step  3a).  If,  in  the  above,  one  has 
x(i,  j)  <  c(i,  j)  ,  node  j  can  be  labeled  if  either  L^(a)  or  L^a)  is 
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of  the  form  b  or  if  c(l)  >  0  (Steps  3b  and  3c).  On  the  other  hand,  if 
L  and  j  are  both  labeled  and  L^(j)  -  b+,  l^Q)  ■  *,  ^(j)  ■  0,  and 
x(J  ,  L)  >  0,  then  L.2^  j)  may  be  assigned  the  value  i  (Step  3d),  thereby 
allowLng  one  to  label  nodes,  k,  from  j  If  x(j  ,  k)  <  c(j  ,  k)  .  These 
poLnts  are  the  basis  of  the  maximum  flow  algorithm  presented  below. 

MAXIMUM  FLOW  ALGORITHM 

1.  Set  c(a)  ■  c(a)  ,  all  a,  and  all  x(l,  j)  «■  0. 

2.  Assign  the  source  the  label  [-,  *,  ®]  and  let  all  other 
nodes  be  unlabeled. 

3.  If  S  is  unlabeied* 

a.  i  is  labeled,  j  unlabeled,  and  x(j ,  i)  >  0,  then  assign 
j  the  label  [i  ,  *,  min  (x(j  ,  l),  L.j(i))]. 

b.  i  is  labeled  with  L^(i)  or  L^(i)  equal  to  b  ,  j  is 
unlabeled,  and  x(i,  j)  <  c(l,  j) ,  assign  j  the  label 
[l+.  *,  min  (c(i,  j)  -  x(l,  j)  ,  L-j(  i) )  ]  • 

c.  I  is  the  source  or  I  is  labeled  with  L^(l)  =  b+  and 
c(l)  >  0,  j  is  unlabeled,  and  x(i,  j)  <  c(l,  j) ,  assign 

j  the  label  [i+,  *,  min  (c(l,  j)  -  x(i,  j)  ,  L3(i)  ,  c(i))]. 

d.  I  is  labeled,  j  is  labeled  with  c(j)  =  0,  L^(j)  =  b+, 

L^O)  “  *»  and  x(j  ,  l)  >  0,  change  the  label  at  node  j 
to  [Lt  (j)  ,  i",  min  (x(j,  i)  ,  L3(i)  ,  L3(j))]. 

4.  Repeat  Step  3  until  either  S  is  labeled  or  no  more  labels 

can  be  assigned.  In  the  former  case,  increase  the  source-sink 
flow  by  F  ■  min  [c(S) ,  L3(S)]  using  the  flow  augmenting 
routine  below.  Otherwise,  go  to  Step  6. 

•k 

For  undirected  networks,  one  considers  arcs  (i,  j)  and  (j,  i)  at 
the  same  time  in  Step  3  by  taking  advantage  of  the  fact  that  at  most  only 
one  of  x(i,  j)  and  x(J ,  i)  need  be  non-zero. 
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5.  If  c(S)  >  G,  go  back  to  Step  2.  Otherwise  go  to  Step  6. 

6.  Terminate  as  the  source-sink  flow  is  optimal. 

When  the  sink  is  labeled,  the  flow  augmenting  routine  traces 
backwards  a  source-sink  chain  and  sends  min  [c(S) ,  L^(S)]  units  of 
flow  along  this  chain.  The  basis  for  the  path  tracing  routine  is 
the  fact  that  if  node  t  is  labeled,  then  the  chain  found  by  the  label 
ing  procedure  to  node  i  is  the  chain  found  to  either  node  | L ^ ( t)  |  or 
|L2(l)|  followed  by  arc  ([^(1)1,  i)  ,  (i,  1^(1)  |),  or  (l,  |L2(i)|) 
as  the  case  may  be.  (If  L^(i)  =  b* ,  then  | L ^ ( i)  |  *  b.)  If  the  last 
arc  in  this  chain  is  (i,  |  L  ^  ( i)  | )  or  (i,  jL2(i)|),  then  one  may  use 
the  chain  to  | ( i)  |  or  | L^C f)  |  designated  by  L^(L^(l))  or  L^L^l)). 
This  is  the  reason  for  setting  ^(a^j)  *  *  in  Cases  2  and  3.  Also, 
in  Case  3,  L^a^)  is  set  equal  to  *,  since  decreasing  flow  along  arc 
(a^,  ajc+^)  brings  the  flow  at  node  a^  below  capacity,  thus  allowing 
one  to  use  the  path  to  node  a^  designated  by  L^(a^)  in  future  calcula 
tions . 


FLOW  AUGMENTING  ROUTINE 


1 .  Let  F  -  min  [c(S),  L^(S)].  Set  k  «  l,  a^  ■  S,  and  reduce 
c(S)  by  F/2  units . 

2.  Case  1:  L^a^)  *  *  and  Ll^ak^  =  b+'  Set  3kfl  "  b’  ^en 

increase  the  value  of  x(a^_^,  a^)  by  F  units  and  reduce  the 
value  of  c(a^)  and  c(a^^)  by  F/2  units. 


Case  2:  L^a^)  *■'  *  and  L^(a^)  =  b  .  Set  a^^  ■  b. 


Then  reduce  the  value  of  x(a^>  a^j)  ty  ^  units,  increase 


the  value  of  c(a^)  and  c(a^^)  by  F/2  units,  and  set 


L2(akfP  “ 
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V.  JUSTIFICATION  OF  THE  FLOW  AUGMENTING  ROUTINE 
This  section  is  to  show  that  upon  termination  of  the  flow  augment¬ 
ing  routine,  the  constraints  of  the  linear  program  of  Sec.  II  hold. 

The  next  section  will  show,  in  addition,  that  the  algorithm  terminates 
with  the  objective  function  v  maximized.  The  justification  of  the 
flow  augmenting  routine  will  be  broken  into  three  parts.  The  first 
will  show  that  the  node  conservation  equations  are  satisfied;  the 
second  that  the  nonnegativity  and  capacity  constraints  on  the  arc  flows 
are  satisfied;  and  the  third  that  the  node  capacity  constraints  are 
satisfied. 


Before  proceeding  to  these  three  parts,  the  following  theorem, 

which  will  be  used  throughout,  must  be  proven. 

Theorem  3.  Let  Lfc  =  L^a^  if  L^a^  =■  a^  ,  or  Lfc  -  lf 

L2(ak)  =  a^^.  Then  L^^  must  have  been  assigned  prior  to  L^. 

Proof.  By  Step  3  of  the  algorithm,  L^(a^^)  must  have  been 

obtained  prior  to  L.  ,  and  the  theorem  holds  if  L.  .  .  ■*  L,  (a,.,). 

r  k’  k+l  1'  kfl 

Suppose  Lj^  ■  ^2^akfP  ’  ^en  c^aicfp  “  0  by  Step  3d.  Furthermore, 
in  determining  a^^  from  the  flow  augmenting  routine,  Case  1  must 
have  occurred,  for  if  not,  ^(a^^)  is  immediately  set  to  *  contradic¬ 
ting  =  L0(a^^).  Thus  L^Ca^)  *  a^p  and  could  only  have  been 

assigned  this  value  by  algorithm  Step  3b,  proving  the  theorem. 


CONSERVATION  CONSTRAINTS  HOLD 

The  node  conservation  constraints  (relationships  III)  are  shown 
to  hold  by:  Theorem  A,  which  shows  that  the  procedure  terminates  with 
no  intermediate  nodes  in  the  path  traced  being  the  source  or  the  sink; 
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Lemma  2,  which  shows  that  flow  may  only  leave  the  source  and  may  only 
enter  the  sink;  and  Theorem  5  which  shows  that  the  conservation 


equations  are  indeed  satisfied. 

Theorem  4.  The  flow  augmenting  routine  terminates  with  a^  =  S. 
Furthermore,  if  a^  a^  is  the  sequence  obtained,  then  k  4  l 

implies  a^  /  S  and  k  ^  n  implies  a^  /  S. 

Proof.  Suppose  a^,  k  >  1,  is  not  the  source.  Since  either  a+ 
is  equal  to  L,(ak  ,),  or  aR  is  equal  to  L  (a^)  or  L^a^)  ,  it 
follows  from  the  steps  of  the  labeling  routine  that  a^  is  labeled  and 
hence  L^(ak)  ,  and  possibly  L^(ak)  ,  exist.  Thus  a^^  exists.  Further¬ 
more,  if  all  a^  ^  S,  the  sequence  of  a^  must  be  infinite,  and  Theorem 
3  would  be  contradicted.  Thus  the  flow  augmenting  routine  terminates 
with  k  equal  to  some  number,  n.  By  Theorem  3,  a^  =  S  implies  k  =  1, 
since  the  algorithm  terminates  once  the  sink  is  labeled;  and  a^  =  S 
implies  k  =  n,  since  the  flow  augmenting  routine  terminates  once  the 
source  is  reached. 

Lemma  2.  x(j ,  S)  =  0  for  all  j  $  B(S)  and  x(S,  j)  *  0  for  all 
j  e  A(S). 

Proof.  Suppose  the  proposition  holds  at  the  beginning  of  the 

flow  augmenting  routine  and  a. . a  is  the  sequence  of  nodes 

1  n 

found.  L  .  =  S+,  where  L  .  is  as  defined  in  Theorem  3,  for  if  not, 
n-1  n-l 

either  or  ^2^an-P  eclua*  t0  ^  »  implying  x(an  S)  >0. 

Thus  flow  is  increased  on  arc  (S,  a  .) ,  and  since  only  a  =  S,  flow 

n- 1  n 

is  unchanged  on  ail  other  arcs  with  S  as  an  endpoint.  Similarly, 

L^  =  a*,  and  flow  is  increased  on  arc  (a^,  S)  and  unchanged  on  all 
other  arcs  with  S  as  an  endpoint.  Thus  it  holds  at  the  routine's 
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sr 


conclusion  and,  consequently,  the  next  time  the  routine  is  entered. 
Since  it  holds  Initially  with  all  x(l,  j)  =  0,  the  theorem  follows 
from  induction. 

Theorem  5.  At  the  conclusion  of  the  flow  augmenting  routine,  the 
following  node  conservation  constraints  (III),  repeated  below,  hold: 

HI  a)  x(S,  j)  x(j,  S)  =  v 

JeA(S)  jcB(S) 

b)  ),  x(i,  j)  -  Y  x(j,  i)  =  0  i  J  S,  S 

JcA(i)  jeB(i) 

c)  £  x(S,  j)  -  ^  x(j,  S)  =  -v. 

jeA(S)  JcB(S) 

Proof:  Suppose  it  holds  initially.  Let  a  .  a  be  the 

sequence  found  by  the  flow  augmenting  routine.  In  Step  2,  x(ak  ak) 
is  increased  by  F  units  if  Case  1  holds,  and  x(ak>  a^)  is  decreased 
by  F  units  if  Case  2  or  Case  3  holds.  In  all  three  cases, 

£x(ak>  j)  _  ak)  »  J  e  A(ak)  ,  j  it  B(ak>  ,  is  decreased  by  F 

units,  and  £x(a ^  ,  j)  -  £x(j\  a^)  ,  j  ,  A(ak+1>  ,  J  e  B(afcf  J  , 

is  increased  by  F  units.  Consider  the  expression  £*(a.  ,  j)  . 

'&(}>  »  J  e  A(a^)  ,  J  e  BCa^  .  If  1  <  i  <  n,  this  expression 

is  increased  by  F  units  when  k  =  i  -  l ,  and  decreased  by  F  units  when 
i ,  f o r  a  total  increase  of  zero.  For  i  =  1,  the  expression  is 
decreased  by  F  units  when  k  =  1;  and  for  i  -  n,  the  expression  is 
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increased  by  F  when  k  =  n  -  1.  Since  a.  ■  S  and  a  =  S,  the  theorem 
J  In 

holds  at  the  conclusion  of  the  routLne  and,  consequently,  holds  ini¬ 
tially  the  next  time  the  routine  is  entered,  with  v  increased  by  F 
units.  Since  it  holds  when  the  routine  is  entered  the  first  time  with 
v  =  0,  the  theorem  follows  from  induction. 

Thus  it  has  been  shown  that  the  balance  equations  always  hold  at 
the  conclusion  of  the  flow  augmenting  routine. 


ARC  FLOW  CONSTRAINTS  HOLD 

Using  Lemmas  3  and  4,  Theorem  6  will  show  that  the  flow  augmenting 

routine  preserves  the  arc  constraints,  0  £  x(i,  j)  £  c(i,  j)  . 

Lemma  3 .  Let  L^  be  as  defined  in  Theorem  3  and  let  L^(a^)  be 

the  value  of  L^Ca^)  immediately  after  the  value  of  L^  is  assigned. 

Then  L^(a^)  is  increasing  in  k  and  all  L^(a^)  ^  F. 

Proof:  Let  L„(a,  .  .)  be  the  value  of  L„(a,  ,  .)  immediately  prior 
-  3  k+l  3  k+l  J  r 

to  assigning  the  value  of  L^.  From  Step  3  of  the  algorithm  the  L^(i) 

cannot  increase.  Thus  it  follows  from  Theorem  3  that  L,(a.  ,  .)  ^  L-(a,  .  ,) 

3  k+l  3  k+r 

Also  from  Step  3,  L^(a^)  ^  L^(a^^).  From  Step  1  of  the  flow  augmenting 
routine  F  s  L^(a^)  ,  proving  the  theorem. 

Lemma  4.  No  arc  appears  more  than  once  in  the  chain  designated 
by  the  sequence  of  nodes  a^,  ...,  a^  and  the  L^  found  by  the  flow 
augmenting  routine. 

Proof:  Suppose  the  proposition  is  false,  and  let  a  =  b  be  the 

m 

last  node  in  the  sequence  belonging  to  an  arc  that  appears  more  than 

once.  Then  there  is  a  k  <  m  such  that  a.  =  b.  Furthermore,  from 

k 

Theorem  3  and  Step  3  of  the  algorithm,  L^  =  L^(b)  and  L^  =  L2(b)  . 
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The  existence  of  L_(b)  Implies  c(b)  =  0,  and  since  L^(b)  could  not 


have  existed  when  L  ,  was  assigned,  L  , 

m-1  m-1 


b  .  Thus  (a  ,  ,  b)  is  a 
m-  i 


repeating  arc.  Furthermore,  a  ,  appears  only  once  in  the  sequence, 

m-  L 

for  if  a  ,  *  a  ,  then  r  <  m-1  would  contradict  Theorem  3  and  r  >  m-l 
m-1  r 

would  contradict  our  assumption  on  a  .  Thus  k  =  m-2  and  L,  =  a+  . , 

r  m  k  m-l 

contradicting  the  fact  that  «  L2(a^). 

Theorem  6.  The  flow  augmenting  routine  terminates  with  all 
0  <:  x(i,  j)  5  c(i,  j)  . 

Proof:  Suppose  it  holds  at  the  start  of  the  routine.  Define 

L^(ak)  as  in  Lemma  3.  From  Step  3  of  the  algorithm,  if  =  a^.  , 

then  s  x(ak,  V,)  ;  tf  Lk  =  then  £  cfa^,  ak) 

-  x(a,  . .  ,  a,).  Since  this  holds  for  all  k,  it  follows  from  Lemmas  3 
k+i  k 

and  A  that  increasing  flow  along  the  chain  found  by  F  units  preserves 
the  conditions  of  the  theorem,  and  it  holds  the  next  time  the  routine 
is  entered.  Since  it  holds  initially  with  all  x(i,  j)  =  0,  the  theorem 
follows  from  induction. 


NODE  CAPACITIES  HOLD 

It  will  be  shown  by  Lemma  5  and  Theorem  7  that  the  node  capacities 
are  not  violated. 

Lemma  5:  These  relationships  hold:  c(S)  =  c(S)  -  £\(S,  j)  , 

j  p  A(S)  ;  c(S)  =  c(S)  -  £]x(k,  S)  ,  k  e  B(S)  ;  and  for  i  ^  S,  S, 

c(i)  =  c(i)  -  ^[£x(i,  j)  +  £x(k,  i)  ]  ,  j  e  A(i),  k  e  B(i)  . 

Proof:  Initially  the  proposition  holds  with  all  x(i,  j)  =  0  and 

all  c(i)  =  c(i) .  Suppose  it  does  hold  at  one  point  and  a,,  ...,  a 

In 

is  the  sequence  of  nodes  that  the  flow  augmenting  routine  generates. 

The  value  of  x(a  ,  a  .)  =  x(S,  a  .)  is  increased  by  F  units  since, 
nn-i  n-I 


-22- 


by  Lemma  2,  the  value  of  x(a  ,  ,  S)  cannot  be  decreased  by  F  unLts. 

n- 1 

Also,  c(S)  is  reduced  by  F/2  units  in  Step  2,  Case  l,  and  by  another 

F/2  units  in  Step  4  for  a  total  of  F  units,  and  the  first  equation  is 

preserved.  Similarly,  the  value  of  x(a2>  =  x(a2>  S)  is  increased 

by  F  units,  and  c(S)  is  reduced  by  F/2  units  in  Step  l  and  by  another 

F/2  units  in  Step  2,  Case  l,  preserving  the  second  equation.  For 

k  /  1,  n,  note  that  a^  S ,  S  by  Theorem  4.  Note  also  from  Step  2  of 

the  flow  augmenting  routine  that  if  either  x(a^,  or  x^ak  L  ’  ak^ 

is  increased  by  F  units,  c(a^)  is  reduced  by  F/2  units;  and  that  if 

either  x(a.  ,  a.  or  x(a,  ,  ,  a,  )  is  decreased  by  F  units,  c(a, )  is 
k  k+1  k-l  k  k 

Increased  by  F/2  units.  Thus  the  third  equation  is  preserved  for  a^. 

Since  the  third  equation  is  unaltered  for  all  nodes  not  in  the  sequence, 

a,,  ....  a  ,  the  lemma  follows  from  induction. 

1  n 

Theorem  7.  At  the  conclusion  of  the  flow  augmenting  routine, 
£x(i,  j)  £  c(i),  j  *  A(i)  for  i  i  S  and£x(k,  S)  5  c(S)  ,  k  e  B(S)  . 

Proof :  From  Lemma  5,  it  is  sufficient  to  show  all  c(i)  s  0.  Let 

a^ ,  ....  a^  be  the  sequence  of  nodes  that  the  flow  augmenting  routine 
generates  and  let  L  (a  )  be  defined  as  in  Lemma  3.  Initially  all 

J  K 

c(i)  ^  0,  since  c(i)  >  0.  Assume  it  holds  at  the  beginning  of  the 
flow  augmenting  routine.  The  only  way  to  reduce  c(a^)  ,  i  /  1,  n,  is 
for  Case  1  to  occur  in  Step  2  for  k  =  i  and  k  =  i-l;  otherwise,  the 
value  of  c(a^)  is  raised  by  F/2  units  when  k  is  equal  to  at  least  one 
of  these  values  and  the  largest  net  decrease  possible  would  be  zero. 


If  Case  1  does  occur  for  k  equal  to  i  and  i-l,  then  c(a^)  is  decreased 

by  F  units.  It  also  follows  that  L,(a.  ,)  =  a+  and  L.(a.)  =  a"!".  .  . 

3  1  L-l  L  1'  L  L+l 
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Thiis  L^a.  p  must  have  been  assigned  its  value  by  Step  3c  of  the 

gorithm  and  c(ap  ^  ^-'3 t _  1 )  an<*>  consequently,  c(a.)  5  F.  No 

c(l)  can  be  decreased  more  than  once,  for  if  so,  one  a  =  a  L  = 

k  m  *  k 

Ll(am)  and  Lm  -  L^a^,  which  would  contradict  Theorem  3.  Thus  the 
theorem  holds  at  the  collusion  of  the  subroutine  and  hence,  If  the 

routine  is  entered  again,  at  the  beginning  the  next  time.  The  theorem 
itself  follows  from  induction. 

It  has  now  been  shown  that  all  network  constraints  hold  after 
each  iteration  of  the  flow  augmenting  routine.  It  now  remains  to 
show  that  the  algorithm  terminates  with  v  maximized. 
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VI.  JUSTIFICATION  OF  THE  ALGORITHM 

This  section  shows  that  the  maximum  flow  algorithm  terminates 
with  v  maximized  and  that  the  maximum  flow  equals  the  minimum  cut. 

Note  that  if  all  flows  and  capacities  are  integers,  the  algorithm 
will  only  yield  Integer  values  for  the  L^(i);  hence  each  iteration  of 
the  flow  augmenting  routine  increases  v  by  at  least  one  unit.  Thus 
if  all  arcs  >nd  nodes  have  finite  capacities,  the  algorithm  must  even¬ 
tually  terminate.  A  similar  argument  holds  for  rational  capacities. 
Lemma  6  and  Theorem  8  will  show  that,  at  termination,  flow  is  maximized 
and  is  equal  to  the  value  of  a  cut  set. 

Lemma  6.  Suppose  the  algorithm  terminates  because  the  sink 
cannot  be  labeled.  Consider  the  cut  set  where  X  consists  of  all 
labeled  nodes,  and  Y  consists  of  ail  unlabeied  nodes  and  also  all 
labeled  nodes,  a,  such  that  there  is  an  unlabeled  node,  b,  with 
x(a,  b)  <  c(a,  b)  .  Then  (i)  x(a,  b)  =  c(a,  b)  ,  if  a  e  X  -  Y  and 

beY  -  X;  (ii)  c(a)  *  0  if  a  e  X  ft  Y;  (iii)  x(a,  b)  =  0  if  a  r  Y  and 

b  e  X. 

Proof:  Suppose  aeX-Y,  b  e  Y  -  X,  and  x(a,  b)  <  c(a,  b)  . 

Then,  since  a  is  labeled  and  b  is  unlabeied,  a  e  Y  also  for  a  contra¬ 
diction.  Thus  (  l)  is  proved.  Suppose  now  that  a  £  X  fl  Y.  Let  b  be 
such  that  b  is  unlabeied  and  x(a,  b)  <  c(a,  b) .  Then  c(a)  =  0;  if 

not,  node  b  could  either  be  labeled  by  Step  3b  or  Step  3c  of  the 

maximum  flow  algorithm,  and  (ii)  is  proved.  Suppose  now  that  a  e  Y, 
b  e  X,  and  x(a,  b)  >  0.  Node  a  i  Y  -  X,  for  if  so,  a  would  be  unlabeied, 
b  labeled,  and  consequently  a  could  be  labeled  by  Step  3a.  Thus 
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a  e  X  Cl  Y  and,  consequently,  c(a)  =  0.  Also,  either  L^(a)  *  c  or 
L^(a)  =  c  for  some  c,  otherwise  node  a's  label  could  be  changed  by 
Step  3d.  Let  d  be  an  uniabeled  node  such  that  x(a,  d)  <  c(a,  d)  . 

Such  a  node  must  exist  since  a  e  X  D  Y.  Then  node  d  could  be  labeled 
by  Step  3b  for  a  contradiction.  Thus  the  theorem  Is  proved. 

Theorem  8.  The  maximum  flow  algorithm  terminates  with  v  maximized. 
Furthermore  the  maximum  flow  is  equal  to  the  minimum  cut  separating  S 
and  S. 

Proof:  Since  by  Theorem  2  the  maximum  flow  cannot  exceed  the 
value  of  any  cut  set,  it  is  sufficient  to  show  that,  at  termination, 
v  is  equal  to  the  value  of  some  cut  set.  Suppose  the  algorithm  termi¬ 
nates  with  c(S)  =  0.  Consider  the  cut  set  where  X  =  {S}  and  Y  *  [N, 
the  set  of  all  nodes}.  The  value  of  this  cut  set  is  c(S)  ,  and  by 
Theorem  1  and  Lemma  2  this  is  equal  to  v.  Suppose  the  algorithm 
terminates  because  the  sink  cannot  be  labeled.  Let  X  and  Y  be  defined 

as  in  Lemma  6.  If  S  e  X  0  Y,  then  v  =  c(S)  by  Lemmas  2,  5,  and  6  (ii). 

If  S  e  X  -  Y,  then  by  Theorem  1  and  Lemma  o  (iii),  v  =  y*x(i,  j)  , 
i  €  X  -  Y,  j  s  Y.  This  sum  may  also  be  expressed  as  v  =  ^  x(i,  j) 

+  £  x(k,  l)  -  J]  x(m,  i)  where  i  e  X  -  Y,  j  E  Y  -  X,  k  e  X,  i  ^  X  fl  Y, 

m  e  Y.  The  first  term  in  this  sum  is  the  total  arc  capacities  of 

cut  (X,  Y)  by  Lemma  6  (i)  ,  and  the  second  term  is  the  total  node 
capacities  of  cut  (X,  Y)  by  Lemma  6  (ii)  and  (iii),  Lemma  5,  and 
Theorem  5.  The  third  term  is  zero  by  Lemma  '  (..ii).  Thus  v  is 
maximized  and  is  equal  to  the  value  of  the  minimum  cut  set. 

This  completes  the  justification  of  the  algorithm. 
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VII.  EXAMPLE 


This  section  uses  the  algorithm  presented  to  find  the  maximum 
flow  from  S  to  S  for  the  undirected  network  of  Fig.  2.  For  purposes 
of  labeling,  the  arcs  are  examined  in  the  order  of  their  listing  in 
Fig.  2.  Among  the  results  listed  for  each  iteration  are  the  final 
labels  obtained,  the  path  found  by  the  flow  augmenting  routine,  the 

resulting  arc  flows,  and  the  final  c(i)  for  the  nodes.  Note  that  in 

the  fourth  iteration  (Fig.  6)  the  chain  found  contains  a  cycle.  In 
the  fifth  iteration,  S  was  not  reached,  terminating  the  procedure. 

The  minimum  cut,  determined  by  the  sets  X  =  [S,  b]  and  Y  =  [a,  b,  c, 

d,  e,  S]  ,  consists  of  arcs  (S ,  a)  and  (S,  d)  ,  and  node  b. 

It  was  mentioned  earlier  that  any  flow  pattern  of  value  v  can 
be  expressed  as  the  sum  of  flows  on  directed  chains  from  source  to 
sink  and  directed  cycles,  such  that  flow  is  non-zero  in  at  most  one 
direction  between  any  pair  of  nodes  and  such  that  the  flows  on  the 
chains  sum  up  to  v.  The  flow  pattern  the  algorithm  finds  can  be 
expressed  as  the  sum  of  directed  chains  only.  It  is  fairly  simple  to 
find  such  chains. 

Note  that  the  iterations  of  the  flow  augmenting  routine  yield 
chains  (not  necessarily  directed)  such  that  the  total  flow  along  all 
of  them  sums  up  to  a  maximum  flow  pattern  if  one  relaxes  the  constraints 
x( i ,  j)  £  c(i ,  j)  to  x(i,  j)  -  x(j ,  i)  ^  c(i ,  j) .  It  is  desired  to 
replace  these  chains  by  directed  chains  such  that  at  most  only  one  of 
x(i,  j)  and  x(j ,  i)  is  non-zero. 


Iterat  ion 


Arc 

Flow 

Node 

Label 

c(a) 

(S,a) 

1 

S 

(-  >  *,  «) 

8 

(a,S) 

1 

a 

1 

(S,d) 

1 

b 

1 

(d,e) 

1 

c 

1 

(e,S) 

1 

d 

(S+,  *,  2) 

9 

(b,d) 

0 

e 

(d+,  *,  2) 

9 

(b,c) 

0 

S 

(e\  *,  1) 

8 

(c,S) 

0 

(S,b) 

0 

(a  ,b) 

0 

(c  ,e) 

0 

Fig.  4  --  Iteration  2;  The  chain  found  is  S,d,e,S 


1 


Arc 

Flow 

(S,a) 

1 

(a,S) 

1 

(S,d) 

2 

(d 

1 

(e,S) 

1 

(b,d) 

1 

(b,c) 

1 

(c,S) 

1 

(S,b) 

0 

(a  ,b) 

0 

(c,e) 

0 

Node 

Label 

c(a) 

S 

( -  ,  * ,  ®) 

7 

a 

1 

b 

(d+,  *,  1) 

0 

c 

(b+,  *,  1) 

0 

d 

(S+,  *,  1) 

8 

e 

(d+,  *,  1) 

9 

S 

(c+,  *,  1) 

7 

Fig.  5  --  Iteration  3:  The  chain  found  is  S,d,b,c,S 


Fig.  7  --  Iteration  5;  The  cut  set  is  determined  by  the  node 
subsets  X  =  js,b|  and  Y  =  ja  ,b  ,c  ,d  ,e  ,S  f.  It 
|  consists  of  node  b  and  arcs  (S,a)  and  (S ,b) 

I  > 

i  ■ 

I 

! 

i 

i 
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Suppose  the  flow  pattern  consists  of  a  unit  of  flow  along  the 

chain,  a,,  a  ,  which  reduces  the  flow  on  arc  (a..,,  a.)  by  one 

I  n  l+1  l 

unit.  Then  there  must  also  be  a  unit  of  flow  along  a  chain,  b^  ,  ..., 

b  with  b.  =  a.,,  and  b,  ,  ,  =  a..  One  may  replace  these  two  chains  by 

m  j  l+1  j+1  l  1  r  J 

the  two  others ,  a,  ,  ...,  a.  ,,  b. . , ,  ...,  b  and  b.  ,  ...,  b.  ,,  a.,., 

I  l-1  j  +  1  m  1  j-1  l+1 

.  .  .  ,  a  ,  thus  reducing  each  of  x(a.  ,  a. , and  x(a. , ,  ,  a.)  by  a  unit, 
n  L  L+1  l+1  l  J 

Repeated  application  of  this  technique  eventually  eliminates  all 
wrong  way  flows  and  thereby  yields  directed  chains.  Similar  reductions 
on  directed  chains  which  use  arcs  (i,  j)  and  (j,  i)  yield  directed  chains 
summing  up  to  a  flow  pattern  where  at  most  one  of  x(i,  j)  and  x(j ,  1) 
are  non-zero  for  all  (i,  j) . 

For  the  example  of  this  section,  the  algorithm  yielded  the  four 
chains: 


l . 

s, 

a, 

S 

2. 

s, 

d, 

e, 

S 

3. 

S, 

d, 

b, 

c, 

4. 

S, 

b, 

d, 

e , 

Eliminating  a  unit  of  flow  on  arcs  (d,  b)  and  (b,  d,  from  chains  3 
and  4  gives: 

1 .  S  ,  a  ,  S 

2 .  S  ,  d ,  e  ,  S 

3.  S,  d,  e,  c,  b,  a,  S 

4 .  S  ,  b ,  c  ,  S . 

Then  a  unit  of  flow  on  arcs  (c,  b)  and  (b,  c)  of  chains  3  and  4 
yields  the  desired  chains: 
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